Transport properties of dense deuterium-tritium plasmas 
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Consistent descriptions of the equation of states, and information about transport coefficients of 
deuterium-tritium mixture are demonstrated through quantum molecular dynamic (QMD) simula- 
tions (up to a density of 600 g/cm 3 and a temperature of 10 4 eV). Diffusion coefficients and viscosity 
are compared with one component plasma model in different regimes from the strong coupled to 
the kinetic one. Electronic and radiative transport coefficients, which are compared with models 
currently used in hydrodynamic simulations of inertial confinement fusion, are evaluated up to 800 
eV. The Lorentz number is also discussed from the highly degenerate to the intermediate region. 

PACS numbers: 51.30.+i, 51.20. +d, 52.65.Yy 



Inertial confinement fusion (ICF) has been a long- 
desired goal, since the implosion of deuterium-tritium 
(DT) capsule at reactor-scale facilities is a visible way 
to generate a virtually unlimited source of energy 0, _2| ■ 
In direct-drive ICF, nominally identical Megajoule-class 
laser beams illuminate the frozen DT capsule. To trig- 
ger the ignition and maximize the thermonuclear energy 
gain, high compression of the DT fuel and a high tem- 
perature of the hot spot should be achieved. In the pro- 
cess of compression, the capsule accelerates and starts to 
converge since the laser beams ablate DT shell's surface. 
The DT shell would experience a wide range of thermo- 
dynamical conditions at densities up to several hundreds 
of gram per cubic centimeters (g/cm 3 ) and temperatures 
from a few to hundreds of electron volts (eV) . 

Understanding and controlling the high pressure be- 
haviors of DT fuel are of crucial interest for the success 
of ignition experiments. The equation of states (EOS) 
of DT fuel are essential for ICF designs, because the low 
adiabat compressibility of the material is dominated by 
EOS [3|]. During the imploding process, impurities from 
ablators or hohlraum are inadvertently mixed into the 
fuel, which strongly influence the burn efficiency so as 
the stabilities of plasmas and the mixing rules. Thus, the 
viscosity and diffusion coefficients of DT mixture are im- 
portant parameters in hydrodynamic modeling treating 
interfaces instabilities |4(, as well as in further investi- 
gating multi-physical effects of high-Z dopants. Finally, 
accurate knowledge of the electronic and radiative trans- 
port coefficients in dense DT plasmas are also important 
for precisely predicting all kinds of hydrodynamical in- 
stabilities grown at the ablation front, at the fuel- ablator 
interface or at the hot spot-fuel interface, and play cen- 
tral roles to describe the evolution of the central hot spot 
against the cold fuel [_J. 

To this end, thermophysical properties of DT mixture 
are determined by combinations of first-principles molec- 
ular dynamics (FPMD) and orbital-free molecular dy- 
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namics (OFMD) @, __] using ABINIT package @. Con- 
sidering DT mixture at high densities, Coulombic pseu- 
dopotential with short cutoff radius (r cut — 0.001 a.u.) 
has been built to avoid overlap between pseudocores [_}. 
As a consequence, large cutoff energy E cut =200 Hartree 
are used. In this letter, we report QMD calculations for 
DT mixtures along the 200 to 600 g/cm 3 isochore with T 
point to represent the Brillouin zone and 128 particles in 
a cubic cell, where FPMD and OFMD are performed up 
to 100 and 10000 eV respectively. Dynamic simulations 
are lasted from 8000 to 200000 steps in the isokinetic en- 
semble. Due to the computational limits of finite temper- 
ature density-functional theory approach, the electronic 
structure is calculated up to 800 eV, and a 4 x 4 x 4 
Monkhorst-Pack fc-point mesh is used. In characterizing 
the plasma states, the ionic coupling and degeneracy pa- 
rameters have been used. The definition of the former 
one is Ta — Z* 2 / (ksTa) (the ratio between the elec- 
trostatic potential and the kinetic energy), and the lat- 
ter one is = T/Tp (the ratio of temperature to Fermi 
temperature). The present studied plasmas undergoes 
from strongly coupled (r^ ~ 60) and degenerated states 
(9 ~ 0.01) to the kinetic states (T u ~ 0.01 and 9 ~ 20). 

Concerning the EOS, the difference between FPMD 
and OFMD lies within 2% accuracy below 20 eV, how- 
ever, as temperature increases, the difference can be 
treated as small as negligible. We have constructed the 
EOS as polynomial expansions (error within 2 %) of 
the density and temperature (P — Aij p l T J ) , which 
can be conveniently used in hydrodynamic simulations 
in ICF or astrophysics (Table. HJ. Figure [T] shows the 
the calculated EOS as a function of temperature at sam- 
pled densities. QMD results are compared with both 
of the ideal and Debye-Huckel models The pres- 

sure of ideal model (Pid) can be viewed as the combina- 
tion of noninteracting classical ions and fermionic elec- 
trons. The self-consistent solution of the Poisson equa- 
tion for a screened charges plasmas leads to the well 
known Debye-Hiickel model, where pressure can be ex- 
plicitly expressed as: Pd = Pid — dpp ■ Here, Ad is the 
Debye length. At a density of 249.45 g/cm 3 , the QMD 
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FIG. 1: (Color online) EOS as a function of temperature at 
249.45 g/cm 3 (Inset is the EOS for p = 400 g/cm 3 ). Previous 
results, including ideal model and Debye-Huckel model, as 
well as PIMC, QLMD, and AAB numerical simulations, are 
also shown for comparison. 

pressures are 9% ~ 40% smaller than those given by 
the ideal model, as temperature lies below 100 eV. How- 
ever, the Debye-Huckel model underestimates the QMD 
results by 60 %, at the lowest temperature (2 eV) con- 
sidered. For very high temperatures, beyond Tp, QMD 
simulation results and the two classical methods merge 
together into the noninteracting classical ions and elec- 
trons (with difference smaller than 1%). Furthermore, 
good agreement in EOS is achieved between the present 
results and other simulation methods, such as, Path In- 
tegral Monte Carlo (PIMC), Quantum Langevin molec- 
ular dynamics (QLMD), and average atom model with 
energy- level broadening (AAB) [llH13f . 

TABLE I: Pressure (kbar) expansion coefficients Aij in terms 
of density (g/cm 3 ) and temperature (eV). 
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Previous QMD simulations in determining the self- 
diffusion coefficients and viscosity were limited at low 
densities (equivalent hydrogen mass density up to 8 
g/cm 3 ) and temperatures (10 eV) (T3 - [l7j . The present 
work substantially extends the study on the transport co- 
efficients into hot dense regime, and comprehensive com- 
parisons with one component plasma (OCP) model are 
described. Here, the self-diffusion coefficient (D) is de- 
rived via particle velocity according to the Green-Kubo 
relation. The viscosity (77) is computed from the auto- 
correlation function of the five off-diagonal components 
of the stress tensor: P xy , P yz , P zx , (P xx - P yy )/2, and 
(P yy — P zz )/2. In order to reduce the length of trajectory, 
empirical functions C[l — exp{— t/r)], where C and r are 



FIG. 2: (Color online) (a) FPMD and OFMD results for the 
diffusion coefficients; (b) Comparison between OFMD simu- 
lations and OCP models for the diffusion coefficients; (c) Vis- 
cosity obtained by FPMD and OFMD; (d) Reduced OFMD 
viscosities are compared with OCP models. For the reduced 
viscosity and diffusion coefficients, the present OFMD results 
are denoted by open squares for the density of 200 g/cm 3 
(black), 300 g/cm 3 (red), 400 g/cm 3 (green), and 600 g/cm 3 
(blue). 



free parameters, have been used to fit the integrals of the 
autocorrelation functions [l7|. The fitting procedure is 
effective in damping the variation and produces reason- 
able approximations to r\ (the computed error lies within 
10%). Due to the fitting procedure and extrapolation to 
infinite time, the total error of 20 % has been estimated. 
The uncertainty in self-diffusion coefficient, where an ad- 
ditional 1/y/N advantage is secured by particle averages, 
lies within 2%. 

Self-diffusion coefficient and viscosity for DT mixture 
calculated by FPMD and OFMD as a function of cou- 
pling parameter Ta at a sampled density of 300 g/cm 3 
are shown in the left panels in Fig. [2j As indicated in the 
figures, FPMD and OFMD results for self-diffusion coef- 
ficients are in generally good agreement (error within 5% 
or better). As for the viscosity, the divergence between 
FPMD and OFMD could reach up to ~ 70% at low tem- 
peratures, however, as temperature increases, they tend 
to merge together. In order to compare these results 
with OCP simulations conveniently, the self-diffusion co- 
efficient D and viscosity rj are reduced to a dimension- 
less form: D* = D/oj p a 2 , and rj* = rj/riiMujpa 2 , where 
Up = (Airrii/M) 1 / 2 Ze is the plasma frequency for ions of 
mass (M = 2.5 amu.). Memory- function analysis of the 
velocity autocorrelation function has been used to ob- 
tain the diffusion coefficient of OCP (D* = 2.95F" 1 - 34 ), 
however, this result is not accurate at T < 4 [18J. Re- 
cently, a more accurate fit (valid at 0.5 < T < 200) to the 
OCP simulations has been provided by Daligault [19|. In 
Fig. H] (b), we compare OFMD data with those OCP 
simulations. The general feature of our simulation result 
agrees with that of Hansen et al. [l8| , but visible diver- 
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gence (3% ~ 50 %) is still observed. The accordance of 
the Daligault's fit with the present results is better at low 
temperatures. As increasing temperature, the Daligault's 
OCP fitting overestimates the OFMD results by ~ 40 % 
at Tu = 0.5. The present results for the reduced viscosity 
are shown in Fig. [2] (d), where results from OCP simu- 
lations are also provided for comparison. By performing 
classical molecular dynamics of OCP, Bastea 20] fit 77* to 
the form: 77* = 0.482r^ 2 +0.629r^°- 878 +0.00188r ij . The 
present simulation results for viscosity indicate a similar 
behavior between dense plasma and normal fluid, that 
is, both of bodily movement of particles and the action 
of interparticle forces contribute to the transport of mo- 
mentum. At large coupling, the viscosity increases with 
decreasing temperature (or adding density), where an 
Arrhenius-type relation (77* = 0.1 x e 0.008r M ^ j s 0Dse rved. 
At the intermediate region, 10 < Ta < 50, contributions 
from the two mechanisms vary with similar magnitude, 
leading to a shallow minimum. OCP simulations indicate 
that the minimum lies at around Ta = 25 19], while our 
FPMD and OFMD results suggest the minimum near 
r,-j = 11 and Ta = 15, respectively. When r,-j < 10, as 
in a gas, the bodily movement of particles is predominant 
and, the viscosity increases with temperature. 

Apart from the EOS and ionic transport coefficients, 
the electronic transport coefficients are derived using 
QMD in the following steps. Ten snapshots are directly 
taken from FPMD equilibrium trajectories up to 100 eV, 
however, beyond this value, the structure are extracted 
from OFMD simulations. Then, the dynamical conduc- 
tivity o~(w) = a\{uj) + i(J2(w) is evaluated through Kubo- 
Greenwood formula as averages of the selected configura- 
tions. The dc conductivity (adc) follows from the static 
limit uj — > of <Ji(u>). In the Chester- Thellung version 
21], the kinetic coefficients La are described as: 



£ y = (-l)W / de*(e)(e- M ) (i+i - 2 >(- 



dfje) 
dt 



), (1) 



with /(e) being the Fermi-Dirac distribution function and 
71 the chemical potential. The electronic thermal conduc- 
tivity K e is given by 



1 JLu 



(2) 



The electronic heat conduction as a function of tem- 
perature along the isochore of 200 g/cm 3 up to 800 eV 
is shown in Fig. [3] (a). Results from theoretical mod- 
els, which are currently used in hydrodynamics simu- 
lations for ICF, as well as previous ab initio simula- 
tions, are provided for comparison. Hubbard model 22] 
is valid for a highly degenerate electron system, where 
the electronic states are treated by using independent 
plane waves. The interactions between nuclei are not 
screened by electrons, and atomic configurations are as- 
sumed to be Debye- or OCP-like. Unfortunately, ther- 
mal conditions in ICF often enters the partially degen- 
erate regime, which clearly goes beyond the hypothesis 
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FIG. 3: (Color online) The present results are compared with 
those of classical models and ab initio simulations, (a) Electri- 
cal thermal conductivity; (b) Lorentz number; (c) Electrical 
conductivity; (d) Radiation thermal conductivity. In panels 
(a) and (c) the labels are combined to present the theoretical 
and QMD results for K e and ado- 



of Hubbard model. Lee- More model [23j, which is not 
quite accurate in the degenerate state, uses different for- 
mulas for the electron collision time in solid, liquid, and 
plasma. The Spitzer model (24[, which exhibits a power 
law on the temperature (T 5 / 2 ), has been dedicated to ki- 
netic plasmas. At high temperatures, our QMD results 
and Lee-More model merged together into the Spitzer 
conductivity. From quantum Boltzmann equation cou- 
pled with Ziman theory, Ichimaru [25[ derived hydrogen 
transport coefficients restricted to the moderate coupling 
r < 2. From Fig. 3, clearly, Ichimaru model results over- 
estimate the present data (up to ~ 30%) at the strong 
coupled regime, and then the difference reduces to 10 
%. Remarkably, our QMD data of the electronic thermal 
conductance for DT mixtures shows reasonable consis- 
tence with recent developed self-consistent average-atom 
model (SCAALP) (2(| and precious QMD simulations for 
hydrogen [j| 0] in a wide range. 

Beyond the thermal conductance, the dynamic conduc- 
tivities are also calculated to examine the Lorenz number 
defined as: 



L = 



h'2 



o dc T 7 e 2 ' 



(3) 



The nature of the screened potential, which is responsi- 
ble for the scattering of the electrons, determines 7. L is 
constant in the degenerate and coupled state (7 is equal 
to 7r 2 /3), reducing to the ideal Sommerfeld number. In 
the nondegenerate case, 7 reaches 1.5966 considering e-e 
collisions. In the intermediate degenerate region, it is dif- 
ficult to deduce the electronic thermal conductivity from 
electrical conductance by Wiedemann-Franz law due to 
the fact that there exists no assumption on the 7 value. 
In Fig. [3] (b) we show the 7 value as a function of degen- 
eracy parameter 9 for different densities. The simulation 
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results indicate that the Lorenz number vibrates around 
the Sommerfeld limit in the degenerate region. As 9 in- 
creases, the departure of Lorenz ratio from the ideal value 
towards the lower limit is observed. 

Results for a^c are presented and compared with those 
obtained by classical models or ab initio simulations in 
Fig. [3] (c). It has been demonstrated that Hubbard 
model recovered the Wiedemann- Franz low in the degen- 
erate limit, and thus the electrical conductance are com- 
puted from the thermal conductance. For the Spitzer 
model, the electrical conductivity is obtained by using 
the Lorenz number at the lower limit of kinetic matter 
including e-e collisions. Our QMD results agree well with 
Ichimaru model [HJ and data by Lambert et al. in a 
wide region between the degenerate state and the kinetic 
state. 

Another useful coefficient, which plays an important 
role in the energy transport formulation of ICF and as- 
trophysics, is the radiative thermal conductivity (K r ). 
As usual, it is connected to the Rosseland mean path 
length (iij = -i^ with kr the Rosseland mean opacity - 
RMO) by the following relation |l|: 



K r = ^a B l R (p,T)T 3 , 



(4) 



where a B is the Stefan-Boltzmann constant. In particu- 
lar, the RMO can be in the form: 



1 

KB. 



a(uj) 



du. 



(5) 



where B'{uS) is the derivative with respect to the tem- 
perature of the normalized Plank's function. a(cu) = 
^-yeri(aj) is the absorption coefficient with n(u) the 
real part of the index of refraction. A comparison for 
K r is made in Fig. [3] (d) between our QMD sim- 
ulations and those given by fully ionized plasma and 



Dyson limit. Zeldovich and Raizer have proven that 
[271 ] the RMO of fully ionized plasma can be obtained 
as n^ eal = 0.0U(Z 3 /A 2 )p(k B T)~ 7 / 2 cm 2 /g. In this re- 
lation, Z is the atomic number, A is the atomic weight, 
and k B T is expressed in kilo-eV. The upper limit of RMO 
(k r < 6.0 x 10 3 Z/(Ak B T) cm 2 /g) has been reported by 
Bernstein and Dyson via applying Schwartz's inequality 
to a particular factoring of the integrand in the defini- 
tion of RMO [HI . Our QMD simulation results indicate 
that K e governs the energy transport process in the fully 
degenerate state. As the matter enters the moderate de- 
generate case, K r rises about four order of magnitude as 
the temperature increases from ~ I0 2 eV to ~ 10 3 eV, 
and becomes dominate when exceeds ~ 0.5. 

In summary, we have determined the EOS and trans- 
port coefficients of DT plasmas within QMD simulations 
in the hot dense regime as reached in future ICF exper- 
iments. The wide range EOS has been built from the 
coupled to the kinetic regime to describe pressures up 
to I0 6 Mbar. A clear chain of simulations in comput- 
ing thermophysical properties of hot dense plasma has 
been demonstrated. The ionic diffusion coefficient and 
viscosity have been simulated and compared with vari- 
ous OCP model simulations. The electronic and radiative 
thermal transport coefficients have also been determined, 
clearly showing different weights when crossing different 
pressure-temperature regimes. The ability to simulate 
these parameters in a self-consistent way shown by our 
results opens a new way to validate classical theoretical 
models currently used in hydrodynamical simulations for 
ICF and astrophysics. 
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